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Recent experiments on the alkali-intercalated iron selenides have raised questions about the symmetry of the 
superconducting phase. Random phase approximation calculations of the leading pairing eigenstate for a tight- 
binding 5-orbital Hubbard-Hund model of AFe2Se2 find that a d-wave (Bi g ) state evolves into an extended s ± 
(Ai g ) state as the system is hole-doped. However, over a range of doping these two states are nearly degenerate. 
Here, we calculate the imaginary part of the magnetic spin susceptibility x"(Qi u) f° r these gaps and discuss 
how the evolution of neutron scattering resonances can distinguish between them. 

PACS numbers: 74.70.-b,74.25.Ha,74.25Jb,74.25.Kc 



I. INTRODUCTION 



Initial experimental reports on the alkali-intercalated iron 
selenide materials of nominal composition AFe2Se2 1 indi- 
cated a number of surprising results 2 which apparently dif- 
ferentiated these materials from their iron chalcogenide and 
iron pnictide superconductor cousins. In contrast to these last 
systems, nearly all of which have hole and electron Fermi sur- 
face pockets present simultaneously, angle-resolved photoe- 
mission (ARPES) studies suggested that there were no hole- 
like Fermi sheets for some dopings. 3 Secondly, the resistiv- 
ity of these samples were all remarkably large, suggesting a 
proximity to a metal-insulator transition. Finally, a phase of 
ordered Fe vacancies supporting a large antiferromagnetic or- 
dered moment arranged in a y/E x y/E block spin pattern was 
reported by elastic neutron scattering. 4 The implications of all 
these anomalous features for superconductivity in these sys- 
tems, which occurs at temperatures as high as 3 IK, is not 
presently clear. Indeed, it is not even established whether 
superconductivity coexists with the vacancy ordered phase, 
arises in Fe vacancy-free regions, or is found in disordered 
vacancy regions of the multiphase samples. 

From the theoretical standpoint, various authors con- 
structed models based on the early ARPES data and inves- 
tigated new physics associated with the absence of the hole 
pockets. Within spin fluctuation theory, the consequences of 
this Fermi surface topology was already discussed by Kuroki 
et al., 5 * who pointed out that this band structure is a 2D 
version of 'Agterberg-Barzykin-Gor'kov" nodeless d-wave 
superconductivity^, since the symmetry-enforced nodal lines 
fall between the Fermi surfaces. This idea was adopted in the 
context of the new systems by several authors. Wang et all 
predicted based on a functional renormalization group (fRG) 
calculation that a d-wave state would be favorable in such a 
situation, with ans ± state a close competitor, and suggested 
that the latter possibility was due to the "marginal" hole band 
just below the Fermi level in their calculations. Maier et al— 
performed similar calculations in the random phase approxi- 
mation (RPA) and also found competing d— and s— wave or- 
der for large electron doping, but only d-wave pairing stabi- 



lized in the absence of hole Fermi pockets 9 . These authors 
also found a strong peak in the dynamical susceptibility not at 
the wave vector Q = (ir, tt) corresponding to the nesting wave 
vector of the two electron pockets in the 1-Fe Brillouin zone, 
but rather close to (n, 0.67r), the vector connecting the clos- 
est flat sides of the rather square electron pockets. Mazir>^ 
and Khodas and Chubukov— have discussed the role of hy- 
bridization and the appearance of an s ± state, in which the gap 
changes sign between the hybridized electron pockets. Other 
authors considered weak-coupling models involving proxim- 
ity to or coexistence with simple magnetic stripes, 12 with the 
y/E x y/E block stated or considered similar electronic struc- 
ture in orbital fluctuation pairing models^ 4 -. Strong-coupling 
models predicting fully gapped s-like states have also been 
proposedi^. 

A final piece of important experimental information was 
obtained when inelastic neutron scattering measurements on 
Rb x Fe2- 2/ Se2 performed by Park et alM reported a reso- 
nance in the superconducting state similar to that observed 
in other Fe-pnictide and Fe-chalcogenide superconductors 18 
except that it was observed not at Q = (jr, 0) (in the 1-Fe 
Brillouin zone) but at Q ~ (7r, tt/2), very close to the value 
predicted by Maier et al s Friemel et aL— then observed a 
dispersion of the resonant mode consistent with band struc- 
ture calculations on RbFe2Se2, within the RPA d-wave spin 
fluctuation picture. This represents strong support for the itin- 
erant nature of the mode, and its ultimate origin in the strong 
scattering between the electron Fermi pockets in the doped 
metallic phase. In these experiments, there was no sign of 
the iron- vacancy ordering; these authors concluded that their 
signal was coming from nonmagnetic, vacancy-disordered or 
free phases, and that an itinerant picture of the phenomenon 
was essentially correct. 

It is important to establish experimentally whether or not a 
d-wave state is realized in this system, since there is accumu- 
lating evidence that all other Fe-based superconductors have 
5-wave (albeit probably s±) symmetry. Observation of a dif- 
ferent symmetry in this case could be consistent with the pre- 
dictions of spin- fluctuation theory, which finds an increasingly 
competitive d-wave instability as hole pockets shrink and 
disappear 9 , or it could represent a genuinely new paradigm 



2 



related to the other unusual features of the alkali-intercalated 
chalcogenides. Recently, Xu et al*22 reported an ARPES mea- 
surement on a KFe2-2ySe2 sample with a small Z-centered 
hole pocket. If a d-wave gap is present, it must by symme- 
try possess nodes on such a pocket. However, the authors 
of Ref. [50 were unable to detect significant gap anisotropy 
on this pocket, and concluded that the data were inconsistent 
with d-wave gap symmetry. On the other hand, synchrotron- 
based ARPES has had a great deal of difficulty observing gap 
anisotropy in these systems, even in situations where other 
probes have provided strong evidence for gap nodes (for a 
discussion, see Ref. [2l|). In addition, we argue below that 
it is difficult to reconcile the inelastic neutron scattering data 
with an s-wave state. We therefore regard the question of the 
symmetry of the gap in these systems as open. 

This paper is concerned with exploring ways in which fur- 
ther neutron scattering experiments on similar systems could 
provide information on the gap structure. We are particularly 
interested in how such experiments can distinguish A\ g s- 
wave and B\ g (d-wave) gap structures in systems where the 
underlying electronic structure consists of two electron sheets 
and small hole pockets which could arise upon hole doping. 
When the AFe2Se2 system is doped so as to move the hole 
bands through the Fermi surface, one expects that as the con- 
ventional Fe-pnictide type Fermi surface, with two electron 
and two hole pockets is recovered, the system will make a 
transition from a B\ g (d-wave) state to an A\ g (,s ± -wave) 
state. However, RPA calculations find that the B\ g pairing 
is surprisingly robust and that a near-degeneracy between the 
A\ g and B\ g channels occurs over a finite doping range where 
the hole pockets first appear. It is in this doping range that we 
will examine how neutron scattering can provide information 
on the gap structure. 



II. MODEL 

In the following, we consider the 5 -orbital Hubbard-Hund 
Hamiltonian 

i,£ i,£'<£ 
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where the interaction parameters U, U', J, J' are given in 
the notation of Kuroki et al^. Here we have used spin ro- 
tational invariant parameters U = 0.92, U' = U/2 and 
J = J' = £7/4. The tight-binding Hamiltonian Ho was fit- 
ted to the full DFT band structure of the parent compound 
KFe2Se2, and the splitting between the two d xz /d yz bands 
and the two d xy bands at the T point was artificially enhanced 
to account for the ARPES results 23 as described in Ref. [H 
The Fermi surfaces and orbital weights | (dz \ vk) | 2 for electron 
dopings of n = 0.15 and 0.05 are shown in Fig. [T] Here, t 
is an orbital index with t G (1,2, 3, 4, 5) corresponding to the 
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FIG. 1: (Color online) Fermi surfaces and orbital weights for (a) 
(n) — 0.15 and (b) (n) = 0. In the latter case the filling is such that 
hole pockets appear around the V point. Colors represent dominant 
orbital content: red (xz), green (yz) and blue (xy). 



Fe-orbitals (d xz ,d yz ,d xy ,d x 2_ y 2,d 3z 2_ r 2) and v and k de- 
note the band and wave vector of the Bloch states. The pair 
scattering processes between the f3 Fermi surfaces of Fig. [T^ 
promote B\ g (d x 2_ y 2) pairing while scattering processes be- 
tween the a and /? Fermi surfaces of Fig. [T]) give rise to 
^lgi^) pairing. Note that the square Fermi surface pock- 
ets found here allow for the possibility of nesting at vectors 
away from (tt, tt) between the two electron f3— pockets, but 
that it is the contribution to the real part of the bare suscepti- 
bility xo due to the opposite signs of the Fermi velocities on 
these two nearly parallel Fermi surface edges which provide 
the dominant enhancement of Fa within the RPA. 



Within the spin fluctuation exchange approach, the pairing 
symmetry for a given set of parameters is determined by the 
eigenvector of the leading eigenvalue of the scattering vertex 
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r(k, k 7 ) in the singlet channel, 

r y (k,k') = Re «lf*(k)a^*(-k) (2) 

x [T /lW4 (k,k' )W = 0)] a£(k>£(-k')- 

Here k and k' are momenta restricted to the electron and hole 
pockets k e Ci and k' e Cj, where z and j correspond 
to either the a or /3 Fermi surfaces, and af,(k) = {di\vk) 
are orbital-band matrix-elements. The vertex function in or- 
bital space Yi 1 i 2 i 3 £ A represent the particle-particle scattering 
of electrons in orbitals t\,i± into £2 , ^3 and in an RPA approx- 
imation is given by: 



(k,k» = 



:U S 



The interaction matrices U s and U c in orbital space are con- 
structed from linear combinations of the interaction parame- 
ters, and their forms are given e.g. in Ref . IIU Here, Xa PA and 



Xo are the spin-fluctuation and the orbital- fluctuation parts 
of the RPA susceptibility, respectively. 

Then, the pairing strength 2 ^ X a for various pairing channels 
a are given as eigenvalues of 



dk'n 



27TO F (k') 



T ij (k,k f )g a (k f )=X a g a (k). (4) 



The eigenf unction g a (k) for the largest eigenvalue determines 
the leading pairing instability and provides an approximate 
form for the superconducting gap A(k) ~ g(k). 

At the doping (n) = 0.15 shown in Fig. [T^, the a hole 
pockets surrounding the T point are absent, and the leading 
pairing instability of Eq. © occurs in the B\ g {d x i_ y i) chan- 
nel. However, as holes are added and (n) decreases, the a 
hole pockets appear as shown in Fig. [TJ) for (n) = 0, and the 
leading pairing instability can occur in the A\ g (s ± ) channel. 

Figure [2] illustrates the momentum dependence of the RPA 
spin susceptibility which drives the pairing along with the mo- 
mentum space structure of the two leading eigenfunctions on 
the Fermi surfaces for four dopings. Here one sees that over a 
range of dopings from 0. 1 to 0.0, these two eigenvalues are re- 
markably close. Thus with doping or possibly pressure there 
can be an evolution from a B\ g (d x 2_ y 2) to an A\ Q (s±) state 
with even the possibility of a d + is stat o 26 i 27 . 

As such an evolution proceeds, one would expect to find 
a change in structure of the neutron scattering resonant re- 
sponse. In the following we calculate the susceptibility in the 
symmetry-broken state as 2 ^ 



2N 



(5) 



k,fiu 



where the 4-momenta are q = (q, uj m ) and k = (k, uj n ). The 
normal and anomalous Green's functions are given as 



G»(k) = tW 9 n+ t^ , F"(fc)= A(k) 



(6) 



Here the matrix elements relating band and orbital space are 

M^ 4 (k,q) =o£(k+q)a£ a *(k)a£* (k+q)<£(k) (7) 

and the quasiparticle energies for band v are given by 
E v (k) = yJQik) + A 2 (k). The neutron scattering cross 
section is proportional to the imaginary part of the spin sus- 
ceptibility in the superconducting state 



E V RPA 



(q,cj) . 



(8) 



^3) The multiorbital RPA spin susceptibility is now given by 



xfZseM^) = {x°M [1 - U a X °(q,u)] X } £i 

(9) 

To understand the type of structure that one can expect 
to see in inelastic neutron scattering experiments for such a 
system in the d x 2_ y 2, s ± and possible intermediate d + is 
states, we have calculated Im xi^ 00 ) from Eqs. ([HO for 
some simple parametrizations of A(k). As seen in Fig. [2 
the Ai g gap exhibits nodes on the electron sheets and has a 
gap which we will parameterize using an extended s-wave 
(xs) (cos 

kx H~ cos ky) form. We will also consider a sim- 
pler isotropic s ± gap. The d-wave gap has the expected 
(cos k x — cos k y ) form except that the gap on the a 2 Fermi 
surface has a phase factor of —1. This is because the spin fluc- 
tuation scattering of d yz pairs from both ol\ and to the 
Fermi surface provides the dominant contribution to the inter- 
action, along with the scattering of d xz pairs to the f3i Fermi 
surface. For a d x 2_ y 2 gap this means that the sign of the gaps 
on the regions of the ql\ and a 2 Fermi surfaces where the d yz 
orbital weight is largest will be opposite to that where the d xz 
orbital weight is largest. Thus the anti-phase appearance of 
the d-wave gaps on the 0L2 hole Fermi surface is simply a re- 
flection of the d xz and d yz orbital weights, which as seen in 
Fig. [T] are out of phase by tt/2 relative with these weights on 
the a 1 Fermi surface. 

With these considerations in mind, we have parameterized 
the gaps as follows: 



Af 
A-(k) 

A^ s (k) 



K (10) 

A zy (cos k x + cos k y ) (11) 

A v (cos k x — cos k y ) (12) 

(A^(k)+tA^(k))/V2 (13) 



x{G^(k + q)G v (k) + F"(-k - q)F v (k)} 



Here the Fermi surface sheet dependent gap amplitudes A u 
are adjusted so that the maximum amplitude is 0.05 on 
each sheet. For the ^ gap, A v on the a sheets is pos- 
itive and on the f3 sheets negative. For the d-wave, A a2 
is negative (out of phase) with respect to the A ai gap. In 
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FIG. 2: (Color online) Left column: ^ = 0) for (top to bottom) dopings n = 0.15, 0.09, 0.05, 0.0 electrons. Next two columns: Leading 
and next leading eigenvector at the same dopings. Red and blue colors indicate positive and negative values, while the symbol size reflects 
the magnitude. For the interaction parameters that we have taken, the d-wave and extended s ± states are nearly degenerate for a range of 
dopings with a change from d to s ± pairing occuring between n = 0.05 and 0.0. 
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the k-space integrations, Eq. [51 the gaps are cut off using 
exp ^(<^(k) — /i) 2 /^o) as k moves away from the Fermi 
surface with = 0.1 eV. 

It has been argued that the d-wave cos k x — cos k y for the 
gap on the f3 Fermi surfaces is fragile and acquires nodes once 
the Se mediated hybridization between the electron pockets 
is taken into accoun t 10 ! 11 . However the nodal regions do not 
make a significant contribution to the spin resonance and the 
simple d-wave form of Eq. (IT2l) provides a suitable represen- 
tation of the Big gap for the 5 -orbital one Fe per unit cell 
model. The hybridization can also alter the usual A\ g s ± gap 
giving rise to an A\ g gap which changes signs between the 
hybridized f3 pockets. The neutron scattering from this state 
is expected to be similar to that of the extended s-wave case, 
Eq. £□]). 

To illustrate how the gap structure is reflected in the inelas- 
tic neutron scattering, in Fig.[3]we show the neutron scattering 
response x"(g,ce;) for an electron doping n = 0.09. This is 
a doping in the range where the A± g and B\ g instabilities are 
close to each other and where small changes can lead to an 
evolution of one state into another or to even the formation 
of an xs + id state. Thus it represents an interesting testing 
ground to explore the neutron scattering resonances. Results 
for momentum transfers q\ = (ir, 0.1 6ir) and q 2 = (ir, 0.7ir), 
corresponding to the peaks in the dynamical sysceptibility of 
Fig. [2 are shown in Figs. [3^ and b, respectively. The s ± 
gapped phase shows the expected resonance response for a 
momentum transfer q\ in which an electron is scattered from 
the a 2 Fermi surface where the gap is positive to the f3i Fermi 
surface where the gap is negative. The extended s-wave xs 
gap shows a similar resonance which arises from scatterings 
between the a 2 and the upper part of the f3i Fermi surface 
where the xs gap is negative. For scatterings with a momen- 
tum transfer q 2 = (7r, 0.697r) shown in Fig. [SJ), the sign of 
the s ± gap is the same on f3\ and f3 2 and the vanishing of 
the coherence factor suppresses the resonance. However, for 
the extended xs-wave, there is a weak response, which for a 
slightly shifted value of momentum transfer q ~ (tt, 0.667r) 
shows a resonance reflecting scattering from the bottom "pos- 
itive" gap region of the f3 2 Fermi surface to the top of the f3i 
Fermi surface where the xs gap is negative. The Brillouin 
zone intensity difference plots which will be discussed below 
provide another way to see this. 

Turning next to the d-wave gap, Fig.[3j) shows a resonance 
for q2 = (ir, 0.697r) associated with the scattering between the 
electron f3i and f3 2 Fermi surfaces. Interestingly, at a lower 
energy there is also a resonance for q 2 = (7r,0.167r). This 
resonance at uo ~ A corresponds to the smaller peak in the 
susceptibility seen for n = 0.09 in Fig. O Here the d-wave 
gap on part of the a 2 hole sheet is out of phase with the gap 
on the /3i sheet, so the coherence factor is non- vanishing and a 
resonance can appear. The xs + id gap also shows a resonance 
response for both qi and q 2 . However, in this case, the energy 
of the resonance at q\ is closer to that of the xs gap and at q 2 
closer to that of the d-wave gap as one would expect. 

Important information helping to identify the gap symme- 
try and structure can be obtained from the qualitative way in 
which the resonance structures at the different wave vectors 





FIG. 3: (Color online) Comparison of neutron response of super- 
conducting states given in Eq. JTol-[T3l) for single doping (n) = 
0.09 (fi = -0.0305). Left: x"(q = O, 0.16tt), lj); Right: 
X"(q=(7r,0.697r),a;). 



disperse with energy. In Fig. [4] we show a set of plots of 
x"{#iw) f° r m e various gap structures given by Eqs. (ITT)1- 
(TT31) along a cut with q x = tt. These susceptibility maps pro- 
vide an alternative way of looking at the response from the 
fixed momentum, variable frequency plots shown in Figs. [3^ 
and b. The map for the s ± gap exhibits structure in the 
(0,7r) and (tt,0) regions associated with scattering from the 
a hole Fermi surfaces to the f3 electron Fermi surfaces. In 
the simple isotropic ^ case, this resonance is near (tt, 0.167r) 
and uj w 1.75A . At higher energies above 2A , a signifi- 
cant broad intensity appears also at larger values of q y due to 
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pair breaking quasiparticle scattering processes between the /? 
sheets. 

The map for the extended xs-wave gap state exhibits a sim- 
ilar resonance peak to the s ± case in the (0, 7r) and (tt, 0) 
regions due to a to f3 scattering. However, as already seen 
in Fig. [3^, the resonance occurs at a lower energy than for 
the s ± gap. In addition, the extended xs state also has reso- 
nances for q ~ (7r,0.667r) and (0.667r,7r). These arise from 
scattering processes between the fa and fa electron sheets, 
whereby an electron scattered from the lower part of the fa 
Fermi surface to the upper part of the fa Fermi surface with 
q no (7r, 0.667r) scatters between regions for which the xs gap 
changes sign. Then at threshold the coherence factor for this 
scattering process is unity and a resonance can appear. In Fig. 
ffl one sees that at higher energies the branch dispersing to- 
wards q y = 7r/2 has larger intensity than the branch dispers- 
ing towards q y = tt. 

The d-wave map shows a strong resonance for uj < 2Ao 
at (tt, 0.697r) and (0.697T, tt) as expected from the uj scan at 
q = (tt, 0.697r) shown in Fig.[3t). This resonance arises from 
fa to fa scatterings, and is seen to have higher intensity to- 
wards q y = 7r as uj increases. In addition, as seen in the uj scan 
at q = (7r, 0.167r), the <i-wave gap also exhibits weaker res- 
onances associated with ol2 to f3 scattering processes. These 
arise because the 7r-phase shift in the sign of the d-wave gap 
on the ot2 Fermi surface gives rise to a resonance for scattering 
an electron from the a 2 Fermi surface to the fa Fermi surface. 

While the g-maps provide a clear way to distinguish be- 
tween an s ± and d-wave state, the extended xs gap and the 
(i-wave gap both show resonances from a-f3 and f3-/3 scatter- 
ing. On the basis of the differences in the spectral weight 
and resonance energies alone, it would be difficult to dis- 
tinguish between xs and d-wave states. However, for the 
xs case, the f3-/3 resonances which onset near (7r,0.667r) 
and (0.667T, tt) show more intensity along the branch dispers- 
ing towards (tt,0) and (0,7r), respectively, while for the d- 
wave case the resonances, which onset near (jr, 0.697r) and 
(0.697r,7r), are stronger towards (tt, 7r). This difference in 
dispersion arises from an interplay between the BCS coher- 
ence factor (1 — ) an d the energy conserving 
^-function S(uj — E^(k + q) — E u (k)), which control the 
scattering phase space. For the xs case, where the resonance 
onsets for q ~ (7r,0.667r), the coherence factor selects scat- 
tering processes in which electrons scatter from the bottom 
of the fa Fermi surface to the top of the fa Fermi surface. 
In this case, as uj increases, the k y separation of the scatter- 
ing states decreases and the xs resonance is stronger towards 
(7r, 0). For the d-wave case, there is a relative change of the 
sign of the gap between the entire fa and fa Fermi surfaces. 
In this case, as uj increases and the k integration for x"((Z> ^) 
runs over a range of different energy cuts of the band energies, 
the peak response shifts towards the commensurate wave vec- 
tor q = (71- , 7r) connecting the centers of fa and fa Fermi sur- 



faces. It is interesting to note that the recent data of Friemel 
et al^ suggest a weak dispersion of the resonance towards 
(7r, 7r), providing additional support for d-wave symmetry. 

Finally, there is the question of identifying the xs + id state. 
While the gap map for an xs + id gap is clearly different from 
the s ± gap and one might differentiate it from the xs case by 
examining the dispersion of the resonances, the difficulty is 
separating it from the d-wave case. That is, while there are 
clearly differences between the xs + id and d-wave gap maps 
shown in Fig.0}), these are quantitative differences which de- 
pend upon the magnitudes of the gaps which are chosen on 
the various Fermi surfaces. The xs + id and d-wave gaps 
can be distinguished qualitatively only at low energies, where 
x"(q, uj) will reflect the full gap in the xs + id case, as op- 
posed to the nodal (i-wave state. For the q-uj maps shown in 
Fig. |4l q x = 7T and the d-wave nodes are not connected by any 
q y . However, nodal excitations can be probed by neutrons at 
other wave vectors. In addition, other experimental probes of 
low-energy quasiparticles can be used to distinguish the d and 
xs + id states. 

While we have presented concrete calculations within a par- 
ticular model and for a given set of interaction parameters, 
leading to resonances at specific energies and wave vectors, 
it is important to stress that the actual resonance energies, 
widths, etc. may differ due to uncertainties in doping, band 
structure, interactions, or other details. In fact, the analysis 
presented here suggests that samples differing in doping by 
very small amounts can have different symmetry order param- 
eters. It is therefore important that experiments like ARPES 
and neutron scattering be performed on the same sample in 
order to draw robust conclusions. Here we have tried to focus 
on differences between various gap symmetries and structures 
which can allow for qualitative distinctions via measurements 
of the neutron resonance at different dopings, and made pre- 
dictions for the evolution of this resonance with doping and 
energy. In conclusion, the AFe2Se2 materials may exhibit an 
Ai g (5 ± -wave), B\ g (d-wave), or possibly xs + id gap. In 
the region where these phases are nearly degenerate, they can 
give rise to resonances associated with a to f3 as well as f3 to f3 
scattering processes. In this case, one can distinguish the A\ g 
state from the B\ g state by examining the dispersion of the f3 
to f3 resonance. 
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